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A\ms. We present a measurement of the dependence of galaxy clustering on galaxy stellar mass at redshift z ~ 0.9, based on the first-epoch data 
from the VVDS-Deep survey. 

Methods. Concentrating on the redshift interval 0.5 < z < 1.2, we measure the projected correlation function, Wp(rp), within mass-selected sub- 
samples covering the range ~ 10' and ~ 10" Mq. We explore and quantify in detail the observational selection biases due to the flux-limited nature 
of the survey, both from the data themselves and using a suite of realistic mock samples constructed by coupling the Millennium Simulation to 
semi-analytic models. We identify the range of masses within which our main conclusions are robust against these effects. Serious incompleteness 
in mass is present below logCM/Mg) = 9.5, with about two thirds of the galaxies in the range 9 < logCM/Mg) < 9.5 that are lost due to their low 
, luminosity and high mass-to-light ratio. However, the sample is expected to be 100% complete in mass above log(M/Mo) = 10. 
^SJ Results. We present the first direct evidence for a dependence of clustering on the galaxy stellar mass at a redshift as high as z ~ 0.85. We quantify 
^ this by fitting the projected function Wp(rp) with a power-law model. The clustering length increases from rg = 2.76*q\1 A"' Mpc for galaxies 
T-H with mass M > 10' Mq to rg = 4.28;^g^j Mpc when only the most massive (M > 10'"^ Mq) are considered. At the same time, we observe a 
I significant increase in the slope, which over the same range of masses, changes from j = 1.67;^jJ jJ^ to y = 2.28;^jj 27. 

' Comparison to the SDSS measurements at z ~ 0.15 shows that the evolution of Wp(rp) is significant for samples of galaxies with M < 10'"'^ M©, 
. while it is negligible for more massive objects. Considering the growth of structure, this implies that the linear bias fo^ of the most massive galaxies 
\^ : evolves more rapidly between these two cosmic epochs. We quantify this effect by computing the value of ht from the SDSS and VVDS clustering 
f"*) amplitudes and find that decreases from 1.5 ± 0.2 at z ~ 0.85 to 1.33 ± 0.03 at z ~ 0.15, for the most massive galaxies, while it remains virtually 
J — constant {b^ ~ 1.3) for the remaining population. Qualitatively, this is the kind of scenario expected for the clustering of dark-matter halos as a 
function of their total mass and redshift. Our result seems therefore to indicate that galaxies with the largest stellar mass today were originally 
• • central objects of the most massive dark-matter halos at earlier times, whose distribution was strongly biased with respect to the overall mass 
^ ' density field. 

Key words. Cosmology: observation - Cosmology: deep redshift surveys - Galaxies: evolution, spectral type - Cosmology: large scale structure 
$H of the universe 



1. Introduction 

In the currently accepted scenario, galaxies are thought to form 
within extended dark matter halos (I White & Reesl[l978l) . which 
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grow through subsequent mergers in a hierarchical fashion. A 
major challenge in testing this general picture is to connect the 
observable properties of galaxies to those of the dark-matter ha- 
los in which they are embedded, as predic ted, e.g., by large n- 
body simulations (e.g. lSpringel et alTl2006h . 



At the current epoch, the observed clustering of galaxies 
found to depend signific antly on t heir s pecific prop- 
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It is somewhat natural to expect that quantities such as lu- 
minosity and, in particular, stellar mass of galaxies should in 
some way be related to the mass of the dark-matter halo. At the 
tim e of observation, thi s is strictly true only for isolated galax- 
ies (ISkibba et al.L |2007|) . but moving back along the evolution- 
ary path of the galaxy one can always find a time when this was 
true, before the galaxy and its halo were accreted by a larger 
dark-matter structure. 

Recent theoretical works seem to indicate that, indeed, there 
exists a fairly direct relationship between global galaxy prop- 
erties (e.g. their stellar or total baryonic mass, or thek lumi- 
nosity) and the halo m ass, before it is accreted by a larger 
dark- matter halo (Minfaii , IConrov et all 120061: IWang et al.L l2006i 
l2007h . In practice, this represents the mass of the dark mat- 
ter halo of the galaxy, at a time when this was still the domi- 
nant central object of its own halo. Sensibly enough, it is Minfaii 
which is defining the global properties of the baryonic compo- 
nent. This complex r elation between galax i es and the (sub)h alos 
in which they reside dSpringel et aliboOUlGao et al.Ll2004 of- 
fers a novel way to interpret the observed clustering properties 
of galaxies. Theoretically, it provides also more direct physical 
insight into the origin of the observed relationships, with respect 
to standard Halo Occu pation Distribution (HOD) models (see 
ICoorav&ShefliLl2002h . 

Observationally, stellar mass has become a quantity mea- 
sured with increasing accuracy, thanks to surveys with multi- 
wavelength photom etry, extending to the near-infrared (e.g. 
iRettura et al.L l2006i) . although some uncertaintie s related to the 
detail ed modelling of stellar evolution remain dPozzetti et al.L 
12007 ^. This makes studies of clustering as a functi on of stellar 
mass possible for large statistical samples. |Li et al.l t2006) have 
measured the dependence of clustering on stellar mass (and other 
properties) in the local Universe, using a sample of ~ 200, 000 
galaxies drawn from the Sloan Digital Sky Survey (SDSS). They 
find, not surprisingly, that galaxies of larger mass are more clus- 
tered than low-mass ones, with the effect increasing above the 
characteristic knee value M* of the Schechter mass function. 
These measurements represent a high-quality reference to which 
high-redshift measurements can be compared to test evolution. 

In studying the clustering of galaxies at different redshifts, a 
selection based on mass should also reduce the typical "progen- 
itor" biases intrinsic when using flux-limited samples to study 
evolution of large-scale structure. Given the strong luminosity 
evolution between reds hift zero and one, and its dependence 
on galaxy spectral type jZucca et all [2006), stellar mass should 
guarantee a more stable parameter on the basis of which com- 
paring galaxies of possibly similar sort. If stellar mass does not 
evolve appreciably below z ~ 1 and the merger rate i s also neg- 
ligible in the same range (which might not be true, see lBell et alj 
( |2006|) ). the growth of clustering observed for objects of fixed 
mass should reflect directly the evolution of the clustering of 
parent dark-matter halos. In fact, the stellar content of galax- 
ies is also thought to evolve between z - I and now, but this 
sh ould correspond to a maximum increase in mass of a factor of 
2 dArnouts et aUl2007l) . A further complication to this idealized 
scenario, comes from the fact that galaxy surveys are inevitably 
flux-limited. Therefore, when extracting mass-selected samples, 
especially at high redshift, particular care has to be taken in con- 
sidering selection effects that might bias the final results, essen- 
tially loosing h igh mass-to-light ratio galaxies at low mass (see 
iLi et alj (l2006l) for discussion and correction of these effects at 
z~0). 

Using the VIMOS-VLT Deep Survey (VVDS), we have the 
possibility to investigate for the first time the evolution with red- 



shift of the clustering in a mass-selected galaxy sample. On one 
side, the VVDS-Deep (see |2.1| below), gives us the opportunity 
to compute accurate stellar masses, by virtue of its multi-band 
coverage that extends to the near-infrared (see iPozzetti et al.l 
(I2007h for detailed discussion on the accuracy and intrinsic limi- 
tations of mass estimates). On the other hand, it provides us with 
an unprecedented area coverage for its depth (Jab ^ 24), thus al- 
lowing the measurement of spatial statistics like the first (mean 
density) and second moments (two-point correlation function) of 
the galaxy distribution. 

Closely related to this study is the work of lPoUo et al.l (|2006|) . 
where the dependence of clustering on galaxy B - band luminos- 
ity at < z >^ 1 has been measured using the VVDS. They show 
that at this epoch, luminous galaxies show not only a higher clus- 
tering than faint objects (similar to what is observed in the lo- 
cal Universe), but that their correlation function is much steeper 
than local counterparts. Since the Z?-band is sensitive to emission 
from young stars, comparison to the results presented in this pa- 
per is of interest to understand how star formation and stellar 
mass assembly relate to the evolution of the local environment 
(and vice-versa). 

In the present paper, we show the first measurement of the 
two point correlation function as a function of stellar mass at 
redshift z ~ 0.85, for galaxies more massive than 10^ Mq. The 
paper is organized as follows. We introduce the data, the way 
stellar masses are derived and our galaxy samples in Sect. |2] 
We present the mock catalogues used throughout this paper in 
Sect. [3] The projected correlation function and its computation is 
introduced in Sect. |4] We discuss the stellar mass completeness 
of our sample and its effect on the measurement of clustering 
properties in Sect. |5] Results are presented in Sect. |6] Summary 
and conclusions are in Sect. [7] 

Throughout this paper we use a Concordance Cosmology 
with Q„, = 0.3 and Qa = 0.7. The Hubble constant is normally 
parametrised via h = i/o/100 to ease comparison with previous 
works. Stellar masses are quoted in unit of /i = 1. All correlation 
length values are quoted in comoving coordinates. 



2. Data 

2.1. The VVDS survey 

The VIMOS-VLT Deep Survey (VVDS) is performed with 
the VIMOS multi-object spectrograph at the ESO Very Large 
Telescope (iLe Fevre et all 120031) and complemented with multi- 
color BVRIJK imaging data obtained at the CFHT and 
NTT telescopes jM cCracke n et all l2003t iLe Fevre et all 12004 
llovino et allboOsTlRadovich et al.l l 



12003'). 



For this work, we use 
the so-called "first epoch" data, collected in the F02 "VVDS- 
Deep" field. This is a purely magnitude limited sample to Iab = 
24, covering an area of 0.49 square degrees with a mean sam- 
pling of ~ 23%. Considering only galaxies with secure (> 80% 
confidence) measurements, this sample includes 6530 redshifts, 
with an rms accuracy (estimated from repeated observations) of 
275 km s"'. This is the same data set used in all previous clus- 
tering analyses of the VVDS. Details about observations, data 
reductio n, redshift measuremen t and quality assessment can be 
found in lLe Fevre et al. I (l2005al) . 

2.2. Mass-limited sub-samples 

Stellar masses for all galaxies in the VVDS catalogue were es- 
timated by fitting their Spectral Energy Distribution (SED), as 
sampled by the VVDS multi-band photometry, with a library of 
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Table 1. Description of the 8 VVDS samples defined in the red- 
shift range z= [0.5- 1.2] 



Sample 


Stellar mass (log(M/Mo)) 


Mean 


Number of 




range 


median 


redshift 


galaxies 


Dl 


9.0- 9.5 


9.25 


0.83 


1221 


D2 


9.5 - 10.0 


9.71 


0.83 


957 


D3 


10.0-10.5 


10.22 


0.84 


670 


D4 


10.5 - 11.0 


10.70 


0.87 


329 


11 


> 9.0 


9.68 


0.84 


3218 


12 


> 9.5 


10.03 


0.84 


1997 


13 


> 10.0 


10.38 


0.85 


1040 


14 


> 10.5 


10.72 


0.87 


370 



O 



I 




-16 



Fig. 1. Comparison of the stellar mass vs. rest-frame B band rela- 
tion in 3 redshift ranges between mock catalogues (black points) 
and VVDS data (red points). 



stellar population models based on iBruzual & Charlo3 (l2003l) . 
Two different histories of star formation were used, a smooth 
one and a complex one with secondary bursts to derive two set of 
galaxy stellar mass. The initial mass function of lChabrieij (l2003b 
was adopted in each case. The resulting typical error on stellar 
mass determination is ~0. 1 dex ( depending on r edshift and stel- 
lar mass). The reader is referred to lPozzetti et al.l ( l2007i) for a full 
description of the methodology used to derive the stellar masses 
and for a discussion of their robustness and intrinsic errors. In 
the next sections, we shall use stellar masses computed using a 
complex star formation history only. In practice, we found that 
the derived clustering properties are not sensitive to the details 
of the method used. 

To measure the dependence of clu stering on stellar mass at 
z ~ 1, similarly to iPollo et al.l (12006 ). we have restricted the 
VVDS-Deep data to the redshift interval [0.5-1 .2]. In this range, 
containing 4285 redshifts, we have constructed a set of mass- 
limited sub-samples, as defined in Tab. [1] These include four 
differentially binned samples, Dl, D2, D3 and D4, with narrow 
mass ranges corresponding respectively to log(M/M0)=[9.O- 
9.5], [9.5-10.0], [10.0-10.5], and [10.5-11.0]. We also consider 
four integrated samples, II, 12, 13 and 14, including galaxies with 
mass larger than a given limit. While the form er are useful fo r 
direct comparison to local SDSS measurements dLi et am2006h . 
results from the integral samples are more robust, given their 
larger size. As we shall discuss in detail in Sect. |5] the lower- 
mass samples are affected, to different degrees, by incomplete- 
ness due to the primary selection in observed flux, and the scat- 
ter in the mass-luminosity relation. A large part of the following 
analysis will be dedicated to understanding this incompleteness 
and its effect on the measured clustering. These tests make ex- 
tensive use of numerical simulations, described in the foUowing 
section. 



3. Simulated mock catalogues 

We have built two sets of 40 mock VVDS surveys, starting 
from semi-analy tic galaxy catalogue s obtained by applying the 
prescriptions of iDe Lucia & BlaizotI (2007) to the dark-matter 
halo merging trees e xtracted from the Millennium simulation 
(Springe l etall l2005h . The Millennium run contains = 2160^ 
particles of mass 8.6 x 10^ /z"' Mq in a cubic box of size 
500/z~' Mpc. The simulation was built with a ACDM cosmo- 
logical model with Q,„ = 0.25, Qa - 0.75, erg = 0.9 and 
Hq - 73km s^^ Mpc~^ . For deta i ls on t he semi-analytic model, 
we refer to iDe Lucia & Blaizot (2007") and reference therein. 
Note that this model uses the Bruzual & Chariot (2003) popula- 
tion synthesis model and a lChabried (120031) Initial Mass Function 
(IMF) to assign luminosities to model galaxies. Our mass es- 



timates are base d on the use of the same population synthesis 
model and IMF dPozzetti et al.Ll2007l) . 

The two sets of 40 1 x1 deg^ hght cone s were generated 
with the code MoMaF dBlaizot et all l2005h for 40 indepen- 
dent lines of sight. The first set is complete in redshift up to 
z ~ 1.7 and contains all galaxies irrespective to their apparent 
magnitude or stellar mass, down to the simulation limit (that 
corresponds roughly to a stellar mass of 10^ Mq). We refer to 
these as the full mocks. They will be used to quantify the stellar 
mass incompleteness and its effect on clustering measurements. 
Nevertheless, these mocks can not be used to built observed 
VVDS-like mock catalogues because do not include galaxies 
above z ~ 1.7. On the other hand, the second set is complete 
in apparent magnitude up to Iab - 24 and was generated inde- 
pendently to the first one. We refer to these mock catalogues as 
the 124 mocks and will use them to quantify statistical errors on 
clustering measurement, using the full spectroscopic VVDS-02h 
observing strategy. 

The Millennium simulation is particularly appropriate as a 
test-bench for the VVDS data, as it is able to reproduce a number 
of basic properties of the data, in particular, the number counts in 
various bands. Additionally, when the VVDS selection function 
is applied, the average redshift distribution of the mock samples 
is in very good agreement with the VVDS. 

Most importantly for the analysis of this paper, the mass- 
luminosity relation of the simulated samples matches fairly well 
that of the data. This is shown in Fig. [1] in three redshift bins 
within the studied range. The general agreement between the en- 
velopes of the observed and simulated relationships is quite en- 
couraging, given the complexity and approximated nature of the 
semi-analytic models. There are clearly some differences, which 
will have to be taken into account when using the simulation 
to interpret some of the effects that are present in the data. The 
most evident one is an apparent excess of massive and lumi- 
nous objects in the mock catalogues, in particular in the most 
distant redshift bin. This will not affect our conclusions when 
using the simulations to understand the level of incompleteness 
in mass due to the flux limit of the survey. This test depends on 
the objects with high mass-to-light ratio, near the lower lumi- 
nosity limit of the data. From the three panels, we can at least 
say that the scatter in the relations for the data and the mocks 
at these luminosities is very similar. We have good reasons to 
hope, therefore, that the results of our tests based on the mocks 
will provide realistic indications on the completeness of the data. 
A more detailed comparison between properties of mock cata- 
logues and VVDS data will be addressed in future papers. 
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4. Estimating the real-space correlation function 

4.1. Methodology 

The two-point correlation function is the simplest estimator used 
to quantify galaxy clustering, being related to the second mo- 
ment of the galaxy distribution, i.e. its variance. In practice, it 
describes the excess probability to observe a pair of galaxies at 
given separatio n, with respect to that of a random distribution 
(lPeeblesl[T980l) . Here we shall use the redshift-space correlation 
function ^(rp,n) . In this case, galaxy separations are split into 
the tangential and radial com ponen ts, and n, a s account for 
redshift-space distortions dPavis & Peeblesill983h . 

The real-space correlation function ^R(r) can be recovered 
by projecting ^(r^, ;7r)along the line-of-sight, as 



w„{r„) 



Jo 



1/2 



\dy 



(1) 



For a power-law correlation function, ^^(r) = {r/roy', this 
integral can be solved analytically and fitted to the observed 
Wp{rp), and find t he best-fitti ng values of the correlation length ro 
and slope y (e.g. lDavis & Peeblesl(ll983h ). In computing Wp(rp), 
the upper integration limit has, in practice, to be chosen finite, 
to include all the real signal without adding extra noise (which 
is dominant above a certain n). In our previous works, we have 
performed extensive experiments (Polio et al., 2005), and found 
that, for the same data, an integration limit of 7T„,ax - 20 h^'Mpc 
provides the best compromise in terms of noise and systematic 
bias of the result. Using 30 or 40 h^'Mpc does not change the re- 
covered Wp(rp) appreciably, while increasing the scatter. Finally, 
fitting Wp(rp) for sep arations < 10 ' M pc using the procedure 
described in detail in iPoUo et all (l2005h . that takes into account 
the covariance matrix of the data and in particular the fact that 
the bins are not fully independent, provides a best-fitting value 
of ro and y for each mass-selected sub-sample. 

To measure in practice ^{ri,,n) from each data sample, we 
use the standard estimator of lLandv&Szalavl(II993h : 

„ , NR{NR-\)GG{rp,n) Nr - I GR(rp,7:) 
' Nc(Nc-l) RR(rp,7T) Ng RR{rp,n) 

where Nq is the mean galaxy density (or, equivalently, the total 
number of objects) in the survey; Nr is the mean density of a cat- 
alogue of random points distributed within the same survey vol- 
ume and with the same selection function as the data; GG{rp,ji) 
is the number of independent galaxy-galaxy pairs with separa- 
tion between and + drp perpendicular to the line-of-sight 
and between n and n + dn along the line of sight; RR(rp, n) is the 
number of independent random-random pairs within the same 
interval of separations and GR{rp,n) represents the number of 
galaxy-random cross pairs. 



4.2. Correction of observational biases 

To properly estimate the correlation function from the VVDS 
data, we need to correct for the different sampling rate, which on 
average is 22.8% but varies with the position on the sky due to 
the VIMOS foot-print and the superpositi on of mult i ple pa sses 
dLe Fevre et all l2005al) . To this end, in Polio et al.1 jlOOSh we 
developed and tested a correction scheme based on computing a 
specific local weight around each galaxy with measured redshift. 
The weight is computed utilizing the angular information of the 
"missed" galaxies contained in the survey parent photometric 
catalogue, which should have the same redshift distribution as 



the spectroscopic sample. In our specific case, in which we are 
considering the clustering of galaxies within a given mass range, 
the parent photometric catalogue should contain only galaxies in 
the same range. This information is a priori not available, as to 
compute it one needs to know galaxy distances. As a remedy to 
this, we used the photometric redshifts that have been measured 
to very good accuracy (crAz/(i+r) = 0.029 for our redshift range), 
over the whole VVDS-02h~Deep field (I lbertet al., 2006), thanks 
to the combination of CFHT-LS and VVDS photometry. This 
uncertainty on the galaxy distance corre sponds to an error o n 
the computed stellar mass of ~0.16dex jPozzetti et al.L l2007l) . 
The chosen redshift range is large enough to assume that edge 
effects due to photometric redshift errors are marginal. We stress 
that here we are not using the photometric redshifts directly into 
the correlation estimate, but only in the correction weight, to 
assign the galaxy to the appropriate (broad) redshift slice and 
mass interval. We used this information to build, for each mass- 
selected spectroscopic data set, a corresponding mass-selected 
parent photometric catalogue containing galaxies in the range 
0.5 < zphot < 1.2. This method allowed us to simplify the 
weight expression ori ginally used for the overall flux-limited 
sample (equation 12, Polio et al.L l2005h . such that the the i-th 
galaxy will have a weight 



np(i)(np(i) - 1) 

\«.(0(«.(0-l) 



(3) 



In this expression. Hp is the total number of galaxies in the 
parent catalogue inside a circle of fi xed radius [choseri to be 40" 
on the basis of mock experiments dPoUo et al.L I2005h 1. centred 
on the i-th galaxy, while n , is the number of spectroscopically- 
observed galaxies inside the same area. The idea behind this 
weighting scheme is that a pair of galaxies will contribute to 
the correlation function with a total weight WjWj, proportional 
to the ratio of the expected-to-measured numbers of pairs within 
the chosen radius. 

Since the sub-samples analyzed in this work are essentially 
volume-limited (above the mass completeness limit discussed in 
the next section), we do not need to apply any other weighting 
scheme as usua lly done for purely flux-limited surveys [as e.g. in 
iLiet al.l(l2006l) 1. 



4.3. Systematic and statistical errors on correlation estimates 

To test for biases in our correction scheme and to quantify statis- 
tical errors on the estimates of Wp(rp), we used the 40 124 mocks 
catalogues described in Sect. [3] From the full 1x1 deg^ light- 
cones, we created 40 fully realistic VVDS-02h surveys, by ap- 
plying the entire selectio n function a n d obse rving strategy of the 
real data as described in iPoUo et all (120051) . For each of the 40 
cones, therefore, we had the possibility to measure both the "true 
Iab ^ 24" and "observed" values of Wp(rp) and the correspond- 
ing power-law-fit parameters to the spatial correlation function. 
Comparing the mean and rms of the difference in the measured 
parameters, we find that (on the average, over the four mass 
ranges), the recovered value of ro is 5% accurate (the maximum 
systematic effect found in the case of the Il-like set of mocks, 
is 10%); similarly, the slope y is consistent with no systematic 
bi as, and a scatter o f 4%. These conclusions are similar to those 
in lPollo et alj ( 12005 ). Finally, the root-mean-square variation of 
Wp(rp) among the 40 "observed" mocks provides us with an es- 
timate of the error bars on our measurements of Wp(rp), ro and y. 
Assuming the simulations are a realistic representation of reality. 
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these error bars also include cosmic variance from fluctuations 
on scales larger than the sample size. 

5. Effects of mass incompleteness 

The VVDS flux limit in the / band translates at different red- 
shifts into different lower luminosity limits. This, in turn, trans- 
lates into a broad mass selection cut at each redshift, reflect- 
ing the scatter in the Mass-Luminosity relation. If not treated 
appropriately, this introduces a bias against high mass-to-light 
(ML) ratio galaxies, i.e. objects that would enter the sample if 
this were purely mass selected, but that are excluded simply be- 
cause they are too faint to fulfill the apparent magnitude limit. 
Furthermore, at a fixed luminosity in the observed band, it will 
be the redder objects (typically more clustered), dominated by 
low-luminosity stars, to have the highest ML ratios and thus to 
be preferentially excluded (see Fig. |4]i. If this is not accounted 
for in some way, it will inevitably affect the estimated clustering 
properties, with respect to a truly complete, mass-selected sam- 
ple. It is therefore necessary to understand in detail the effective 
completeness level in stellar mass of the samples that we have 
defined for our analysis. We have explored this issue using two 
complementary and independent approaches based respectively 
on the data themselves and on the mock samples. The latter will 
also allow us to test directly the effects of the incompleteness on 
the measured clustering; the results of this test will be presented 
in Sect. 



5.1. Observed scatter in the mass-luminosity relation 

A first way to understand the incompleteness in mass of our sam- 
ples is to measure directly the scatter in the mass-luminosity re- 
lation near the flux limit of the survey and from this estimate 
the fraction of missing objects, at the given mass limit. Since 
the flux limit of the survey translates into a diff'erent luminosity 
limit at different redshifts, we need to do this in narrow redshift 
ranges as to track the change of the incompleteness within the 
broad redshift range of the sample (and also take into account 
any potential evolution in the Mass-Luminosity relation). The 
method is illustrated by Fig. |2] where we plot the value of B- 
band luminosity vs. the computed stellar mass for all galaxies in 
the redshift rangez=[0. 8-0. 9]. The "completeness" limit, defined 
as the value of stellar mass above which virtually all masses, 
given the observed distribution, are included in the sample is 
estimated as follows: 1) given a redshift bin (Az = 0.1), we com- 
pute the minimum galaxy luminosity detected at that distance, 
defined as the value of absolute magnitude above which we have 
95% of the galaxies, Mmm, this is given by the long-dashed hor- 
izontal line of Fig. |2] 2) We then consider a luminosity interval 
[M,„i„; M„,ii, - 0.35], chosen as the minimum thickness to pro- 
vide us with statistics sufficient to define the properties - mean 
and scatter - of the (logarithmic) mass distribution, as indicated 
by the short-dashed horizontal line in the same figure. We shall 
define for a specific redshift the completeness limit to be: 



Yciz) = {Y(z)} + crriz) , 



(4) 



where Y(z) - log M(z). Note that the completeness limit defines, 
as a function of z, the value of mass above which, given the 
selection properties of the survey, a sample will be better than 
84% complete in mass (corresponding to one-sided Icr limit). 
Note that here we do not have enough data to perform a re- 
fined Gaussian reconstruction as done by Li et al.l d2006) for the 
SDSS. The completeness locus as a function of redshift is plotted 




<: 



00 
o 



-20 



-18 

I 



-16 



0.8 0.85 

0.85 

z 

0.8 



+ 

CO I CO 
CD CD 

I I I I I I 1 . 1 ... I ,.l.-J,, 



_L 



AM„=0.35 



„ l „ LI.. ■!■ . 1 - 1 I I 



z = [0.6-0.9] 



8 9 10 11 

log(Mj [h-= MJ 



Fig. 2. Construction of the completeness limit from VVDS data 
in the narrow redshift range z= [0.8-0. 9]. The central panel shows 
the value of B-band luminosity vs. the computed stellar mass for 
all galaxies in that redshift range. The long-dashed horizontal 
line defines the minimum galaxy luminosity above which stand 
95% of the galaxies. The top horizontal dashed line defines the 
upper limit in which we compute the mean and standard de- 
viation of the stellar mass distribution of galaxies (top panel) 
standing in a luminosity interval of AMb = 0.35. The right-hand 
dashed vertical line defines the 84% completeness Umit in stellar 
mass. 



as the bold (red) solid line in Fig. |3]in the plane redshift-mass. 
Figure [3] indicates that any sample limited to masses larger than 
IO'^-^Mq will be fairly complete in stellar mass over the whole 
redshift range we plan to explore. This completeness l imit is 
consistent with the one derived by iPozzetti et al.l (l2007h in the 
redshift range [0.9-1.2] in the determination of the galaxy stellar 
mass function from the VVDS data (see their Fig. 9). 

5.2. Mass incompleteness from the mock catalogues 

Encouraged by the realistic appearance of the Mass-Luminosity 
relation for our mock samples (Fig.lTJ, we can try and use them 
to independently test the consistency of the incompleteness level 
estimated in the previous section. 

As a first test, we can apply to the mock samples the direct 
technique used on the data in the previous section. The results of 
this, averaged over the 40 mock samples, are shown in terms of 
a completeness curve in the plane z - M in Fig. [3] (dashed lines). 
Interestingly, we see that the 84% completeness level obtained 
from the simulations is very close to that from the data (upper 
curve), although the mean values of the two mass-luminosity re- 
lationships are slightly different. 

Then, let us exploit the full mock catalogues defined in 
Sect. [3j which are deeper than the magnitude limit of the VVDS 
and include galaxies with stellar mass down to 10"^ Mq. It 
is therefore possible to quantify independently the fraction of 
galaxies we miss when applying the VVDS selection function. 
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0.8 0.9 
redshift z 

Fig. 3. Distribution of stellar masses as a function of redshift, 
in the range [0.5-1.2]. The top (bold) solid line, computed from 
the observed scatter in the luminosity-mass relation, defines the 
limit above which samples are better than 84% complete. The 
bottom line defines the 50% completeness limit. The two dashed 
lines represent the same limits but derived from the Millennium 
simulation. 

Table 2. Mean fraction of galaxies in the redshift range z=[0.5- 
1.2] for which 17.5 < Iab ^ 24 among 40 mocks catalogues (se- 
rie SI) which are complete in stellar mass down to 10^ Mq, built 
from the Millenniuni simu lation with semi-analytical model of 
iDe Lucia & BlaizotI (l2007h 



log(M/Mo) 


C 


9.0- 9.5 


0.37 


9.5 - 10.0 


0.70 


10.0 - 10.5 


0.93 


10.5-11.0 


0.99 


> 9.0 


0.56 


> 9.5 


0.79 


> 10.0 


0.94 


> 10.5 


1.00 



provided that the Mass-Luminosity relation and its scatter are 
comparable in the data and in the simulations. We have seen 
from Sect. [3]that this is quite a reasonable assumption. 

We define, for a given mass selection, the completeness C as 
the fraction of objects recovered when the apparent magnitude 
selection 17.5 < Iab < 24 is applied, with respect to the pure 
mass-selected sample. The average values of C among the 40 
full mocks catalogues are reported in Tab.|2]for z= [0.5- 1.2]. 

These values indicate that for 9 < log(M/MQ) < 9.5, our 
sample is strongly incomplete, with an expected value of lost 
galaxies larger than 60%. This represents an average value over 
the z= [0.5- 1.2] redshift range (while we have seen from the pre- 
vious analysis that the incompleteness is actually a function of 
redshift). More encouragingly, the tests also show that an inte- 
gral sample with log(M/M0) > 9.5 is already ~ 80% complete 
and that samples above log(M/MQ) = 10 are better than 93% 
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Fig. 4. Mean luminosity and color distributions, with their 
Icr variance, for galaxies with stellar masses in the range 
log(M/M0) = [9 - 9.5] in the full millennium mock cata- 
logues, for 3 consecutive redshift ranges between z = 0.5 and 
z - 1.2. Solid red distributions correspond to galaxies fainter 
than Iab = 24, while hatched blue distributions correspond to 
galaxies brighter than this threshold. The distribution have been 
normalized to the total number of galaxies in each redshift range. 
The galaxies missed due to the Iab - 24 flux limit function are 
clearly intrinsically fainter and relatively redder. 



complete both in the differential and integral cases. It is interest- 
ing to use the simulated sample to explore which are the prop- 
erties of the "missed" galaxies. Figure |4] shows the luminosity 
and rest-frame (B-I) color distributions of galaxies fainter than 
Iab = 24 but with masses in the range log(M/M0) = [9.0 - 9.5], 
compared to those brighter than this limit (and thus retained 
within the sample). Interestingly, the population missed by the 
VVDS selection function - at least for the mock samples - is 
fairly red galaxies, as one would expect for higher-than-average 
mass-luminosity ratio objects dominated by old stars. 

Finally, considering the differential 84% completeness limit 
in Fig. |3] we note that - once averaged over the redshift range 
- it roughly corresponds to a mass limit log(M/M0) > [9.5 - 
10]. This is consistent with the integral completeness fraction C 
derived in this section using the simulations, for the same range 
of masses (Table|2|i. 

5.3. Effects on the measured clustering 

We can use the 40 full mock catalogues to go one step further 
and directly quantify the effects of the mass incompleteness on 
the clustering measurements, when compared to purely mass- 
selected samples. 

Figure |5] shows the average and scatter over the 40 mocks of 
the ratio of the projected functions measured respectively with 
and without the VVDS flux limit, within the usual mass ranges. 
As expected, the values of Wp(rp) are significantly affected by the 
mass incompleteness, in particular at very small separations. For 
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samples containing masses below logCM/M©) ~ 10, not only the 
amplitude is underestimated by 15% or worse, but also the shape 
of Wp{rp) is affected, with the average ratio becoming smaller 
and smaller below 1 /z ' Mpc. This can be translated in terms of 
a systematic bias in the measurements of ro and y. The results 
of averaging the relative difference Sxjx - {x,rue - Xobs)/^ over 
the 40 mocks, are reported in Tab. [3] in terms of a relative bias 
e{x) =< dx/x >4o on each parameter In all cases the mass in- 
completeness produces a systematic underestimate of both the 
amplitude and slope of ^(r). 

In particular, we have a systematic ~ 25% underesti- 
mate of the correlation length ro in the lowest-mass sample 
(logCM/Mo) = [9.0 - 9.5]). The slope y is less affected (7%); 
this is related to the fact that the fit is performed out to rp - 
10 h"' Mpc, and thus is not so sensitive to the small-scale flatten- 
ing that is evident in Fig.|5] We also see that the systematic effect 
is significantly less severe for the integral samples (e.g. only 6% 
on ro and 3% on y for log(M/MQ) > 9.5). 

On the basis of these extensive tests, we shall consider as 
reliable in the analysis of the VVDS data only results for sam- 
ples with a minimum mass limit of log(M/Mo) = 9.5, with 
some caution also in interpreting the results from the sample D2 
(9.5 < log(M/Mo) < 10), which is still affected by some residual 
incompleteness. 

In the next section, we shall present the measured values in 
general uncorrected for these effects, but will discuss and show 
how they would change assuming a correction for the incom- 
pleteness as derived in Tab. [3] All plotted error bars on the mea- 
sured points will correspond to the statistical uncertainty only, 
as estimated from the variance in the 40 124 mock catalogues. 



Table 3. Average systematic underestimate of the clustering 
length ro and the slope y of the correlation function introduced 
by the mass incompleteness. The values are computed for each 
parameter as e(x) = < 5x / x >4o, where 5(jc) = Xirue-Xobsemd 
and the average is performed over 40 mock samples. 



log(M/Mo) average underestimate 





f(ro) 


e(r) 


9.0- 9.5 


25% 


7% 


9.5 - 10.0 


11% 


3% 


10.0-10.5 


3% 


1% 


10.5-11.0 


0% 


0% 


> 9.0 


13% 


5% 


> 9.5 


6% 


3% 


> 10.0 


2% 


1% 


> 10.5 


0% 


0% 



6. Results 

6.1. Wp{rp) as a function of mass atz~ 0.85 

Figure |6| shows the projected correlation function Wp(rp) com- 
puted for the four samples II to 14 (left panel) and its best fit 
parameters together with their 1-, 2- and 3cr error contours de- 
rived from the 40 124 observed mock catalogues (right panel). 
The measured values of ro and y for each sample are reported 
for convenience also in Tab. |4] including those for the sample 
Dl for which we know that the incompleteness is severe. No 
correction for incompleteness is applied to the plotted values 

Even accounting for the effect of incompleteness and the size 
of the error contours, we observe a clear increase of the ampli- 
tude and the slope of the correlation function for increasing me- 
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Table 4. Clustering length and slope of the correlation function 
Wp{rp) for the different samples in the redshift range z=[0.5-1.2]. 
Values are derived for 0.2 < rp < 21 h^\Mpc. Associated errors 
have been extrapolated from the variance measured among 40 
mocks catalogues for each stellar mass range. Values have not 
been corrected by underestimate due to stellar mass incomplete- 
ness. 
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Fig. 7. Correlation length and slope measured from V VDS data 
at z ~ 0.85 as a function of the median stellar mass of each 
samples. Left and right panels refer respectively to the differen- 
tial and /nfe^rafeii samples. Filed symbols indicate measurement 
non-corrected by stellar mass incompleteness, while open sym- 
bols are. 



dian stellar mass, even if error contours are quite large in the 
highest mass range where the statistics is poor (Tab. [T]i. This 
trend is also quite clear in Fig. [7] where instead we have used 
the simulation results to correct (open symbols) the observed 
values (filled symbols), dividing them by (1 - e(x)) (Tab. [3]). 
Note how especially for the integral samples (right panels), the 
trend with increasing mass is quite robust and significant. We 
thus conclude that already at z ~ 0.9 higher mass galaxies are 
significantly more clustered than lower mass objects, with a sys- 
tematic trend both in the amplitude and slope of the spatial cor- 
relation function. This dependence on stellar mass is coherent in 
general terms with the trend found as a function of luminosity 
(iPolloet al.Ll2006l) . 



6.2. Comparison to SDSS 

This segregation in term of stellar mass at reds hift z ~ 0.85 i s 
also observed at low redshift in the SDSS data dLi et al.Ll2006l) . 
These authors present their result directly in terms of the mea- 
sured amplitude of Wp(rp) on different physical scales because 
of the departure of the correlation function from a pure power 
law observe d by n umerous authors (e.g. iGuzzo et al.L Il99ll; 
IZehavi et al.L |2004|) . In Fig. |8] we plot our measur e ments of 
Wp(rp) at z ~ 0.85 together with those by iLi et all (|2006|) at 
z ~ 0.15 within the three mass ranges in which we have es- 
tablished that the incompleteness does not strongly affect the 
measurements. A reference power-law Wp{rp), corresponding to 
ro = 5 h"' Mpc and y = 1.8 is over-plotted for reference. 
Blue open circles correspond to correcting the observed values 
of Wp{rp) according to the results of Fig.|5] 

Figure|8]not only shows the dependence of galaxy clustering 
on stellar mass at a given redshift but also, at the same time, the 
evolution of the amplitude and shape of the correlation function 
Wp(rp) with redshift for a given stellar mass range. The observed 
evolution is apparently faster for low-mass objects than for mas- 
sive ones: the amplitude of Wp(rp) increases by a factor ~2-3. 
Conversely, in the high-mass range, log(M/M0) = [10.5 - 1 1.0], 
the amplitudes at z ~ 0.85 and z ~ 0.15 are very similar, within 
the error bars. These results hold even when considering correc- 
tions to the projected function accounting for the stellar mass 
incompleteness (blue empty circles, Fig.[8]|. 



6.3. Comparison to model predictions 

So far, we have used the mock VVDS surveys built from the 
Millennium Simulation to assess where our results are reliable 
against the survey selection effects and established the range 
of stellar masses where this is plausibly the case. In this sec- 
tion we would like to exploit more fully their intrinsic sci- 
entific content. As discussed previously, these mock samples 
were constructed by populating dark matter halos from the 
Millennium Run (^ Springel et al.L I20p5|). using th e latest version 
of the Munich semi-analytic model ( De Lucia & Blaizot, 200?)). 
We have shown this model is in quite agreement with a number 
of observational results. As we already discussed, the "ideal" 
light cones were "observed" to produce realistic replicas of the 
VVDS-Deep survey, with the same sampling, field pattern and 
incompleteness. They are therefore ideal for providing a direct 
comparison, under the same conditions, to the Wp(r;,)measured 
from the VVDS data. This is shown in Fig. |9] where the data 
points are compared to the average of the measurements from 
the 40 observed mock samples (solid line), including a +lcr er- 
ror corridor The agreement of the data and the model predictions 
for the redshift range considered, 0.5 < z < 1.2, is rather good 
for all three mass ranges explored. The VVDS points are consis- 
tent with the predictions of the semi-analytic model within Icr 
or better A more comprehensive comparison of these models 
with the clustering of VVDS galaxies as a function of different 
properties (luminosity, color, spectral type), will be presented in 
a separate paper. 



Finally, we also note that a similar dependence on stellar 
mass is also p r edicte d in the hydrodynamical simulations of 
Wei nberg et al.l (|2004|) . in which the steepening of the correla- 
tion function is particularly marked on scales below 0.5/; 'Mpc. 
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Fig. 6. (Left) Measurements of the projected correlation function Wp(rp) of galaxies with different stellar masses: log(M / Mq) > 9.0 
(open blue squares), > 9.5 (filled red squares), > 10.0 (open green triangles) and > 10.5 (filled magenta triangles) from VVDS 
data in the redshift range [0.5-1.2]. (Right) The best-fit parameters (ro and y) with their associated 1-, 2- and 3cr error contours, 
derived from the variance among 40 mock catalogues. In this plot, no correction for the effect of incompleteness in mass is applied 
(see text). Note however that the correction would affect both ro and j and, for example, move the open square (the most affected 
one) in the right panel without modifying the observed trend. 



6.4. Evolution on galaxy bias 

Our results suggest that the apparent evolution of clustering from 
z ~ 1 to the current epoch has been stronger for the less massive 
galaxies than for the most massive ones. In other words, the rel- 
ative bias of massive to less-massive galaxies was larger at z ~ 1 
than it is today. This behavior is also fairly well reproduced by 
the Millennium mocks. 

In the currently standard scenario of galaxy forma- 
tion, galaxies form and evolve inside dark matter haloes 
d White & Reesl 1 19781) . with the first objects appearing inside 
the most massive of these. These are expected to correspond 
to the highest, more rare fluctuations of the overall matter den- 
sity field. On simple statistical grounds, it can be shown that 
these peaks will be more clustered than the gen e ral da rk mat- 
ter distribution, i.e. they will be biased dKaiserl |1984|) . Thus, 
from very simple considerations, if a larger stellar mass corre- 
sponds to a larger hosting dark matter halo mass, then the most 
massive galaxies at high redshift will correspond to the high- 
est peaks and will be highly biased. Strong clustering and high 
bias are observed both for massive spheroidals at redshifts in the 
range z ~ 1.5 - 2.5, selected on the ba sis of their red optical- 
infrared colors (e.g. iDaddi et akl l2003h . and for strongly star- 
forming gala xies at z ~ 3 selec ted from their strong UV Lyman 
break (e.g. .Steidel et aUl2004 . 

Another important piece of information is that, when se- 
lected on the basis of a fixed rest-frame luminosity, galaxies 
show an increase of their average bias as a function of redshift. 
Using the VVDS data, we have in particular shown that the bias 
of galaxies with M/, < -20 + 5 log h increases by a factor ~ 1 .5 



between z - and z ~ 1.5 jMarinoni et al.Ll2005h . However, we 
also know that the typical luminosity of the general population 
of galaxies increases by more than a magnitude over the same 
redshift range, which implies that at a fixed absolute luminos- 
ity we are in fact looking at different populations of objects at 
different redshifts. Given that there is clear evidence for at most 
a fact or of 2 incre ase in the stellar mass between z = and 
z = 1 jArnouts et al ., 2007), the analysis presented here has the 
virtue of addressing the dependence of clustering with respect 
to a more physical and stable property of the galaxy population. 
Even by allowing for a factor of 2 increase in the mass (i.e. as- 
suming that the progenitors of objects with mass M at z ~ are 
galaxies with mass M/2, at z ~ 1), our results would not change 
significantly. 

It is therefore of interest to interpret in a more quantita- 
tive way our findings, by estimating the evolution of the lin- 
ear bias for galaxies w ith similar mass at different redshifts. In 
iMarinoni et al.l (|2005|) we showed that the bias is in fact non- 
linear, at about the 10% level. This has very little efi'ect on the 
present discussion, and we can safely assume here a linear bias 
model for bi, where the galaxy rms fluctuations cr«_j on a given 
scale R can be expressed as 

crR,g = ^lO"r,,„ . (5) 

with cTg ,„ the mass rms fluctuations on the same scale. We also 
make the further assumption that the bias bi is independent of 
scale, which is very reasonable on sufficiently large scales. It is 
customary to express the rms fluctuation of galaxies and mass 
in spheres of radius /? = 8 h^'Mpc, with the value of cr^ ,,, 
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Fig. 8. Comparison of the projected correlation function Wp{rp) from SDSS data at z~0.15 (black cross: lEiet ali \2d0^} to VVDS 
measurement at z~0.85 (magenta filled dots and solid lines) in 3 differential ranges of stellar masses. The green dashed line is a 
power-law reference line drawn with ro = 5 h ' .Mpc and j - 1.8. Error bars on VVDS measurement have been estimated from the 
variance among 40 mock catalogues. Blue open circles indicate VVDS measurements corrected for stellar mass incompleteness. 



1000 E- 



— I I M MM I I M MM n_ I I M MMI I I M MM □ = 

E 9.5^1og(M/Mo)^10 EE 10^1og(M/Mo)^l'0. See 



100 E- 



10 E- 



I I MM I I M MM 

;iog(M/MQ)^ll 




0.1 



1 

[h- 



10 0.1 
Mpc] r 



1 

[h-^ 



10 0.1 
Mpc] r 



1 

[h-^ 



10 
Mpc] 



Fig. 9. Comparison of the VVDS measurements of the projected coiTelation function Wp(rp) as a function of stellar mass (filled 
circles), with that obtained from 40 mock VVDS surveys constructed from the Millennium Simulation. The dashed area describes 
the rms among the 40 simulated catalogues. 



adopted as standard expression for the normalization of the mass 
power spectru m at the present epo ch. Here we adopt the value 
o"8,m - 0.76 (Spergel et al.L l2007l) . This value at z = can be 
scaled to the epochs of interest here by considering the growth of 
fluctuations in our assumed model, cr8_„,(z) = o-^,m(z = O)D(z), 
where D(z) is the linear growth factor of density fluctuations, 
D(z) - g(z)/[g(Q)(l + z)], and the normalized^growth factor g(z) 
can be approximated as jCarroU et al.lll992HMo & Whitell2002l) 



with and Qa the density parameters in non relativistic matter 
and cosmological constant. Note that in this expression. 



nm,0(l +Zf 



SI A = ^a(z) = 



Q 



A,0 



EHz) ' '''' EHz)' 

with Q.m.0 and Qa o their present-day values and 

e\z) = [flA,,) + (1 - f2o)(l + zf + f2™,o(l + zf] , 



(7) 



(8) 



the normalized expansion factor. 
5 I The corresponding value of the rms galaxy fluctuations on 

^(z)^ ^i^m[i^;"-OA + (l+Om/2)(l+QA/70)J , (6) the same scale can be estimated from the parameters of the cor- 
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Table 5. Measured values of the linear bias for the VVDS and 
SDSS samples, within 3 stellar mass ranges (see text). 



Sample 


log(M/Mo) 




b(z) 






SDSS (z 


-0.15) VVDS (z~ 0.85) 


D2 


9.5 - 10.0 


1.26 ± 


0.03 1.29 ±0.10 


D3 


10.0 - 10.5 


1.37 ± 


0.03 1.32 + 0.12 


D4 


10.5 - 11.0 


1.40 ± 


0.03 1.62 + 0.18 



relation function under the assumption of a power-law form, 
^(r) = {r^lry (Peebles, 1980j, as 

where 



' (3-r)(4-r)(6-r)2r 

To be consistent with our work at z = 0.85, we estimated the 
parameters ro and y for the SDSS at z ~ 0. 15 from the published 
measurements of Wp{rp) by Li et al.l (l2006h . fitting over the same 
range used for the VVDS (0.18 < < MSIr^Mpc). The re- 
sulting values of bi (with Q^.o - 0.3 and Qa,o - 0.7) over the 
three mass ranges and at the two redshifts considered are plotted 
in Fig. [To] The values are also reported for reference in Tab.|5] 

Figure [TOl shows that for galaxies in the highest-mass range 
considered here (10.5 < log(M/Mo) < 11), the evolution of the 
linear bias between redshift 0.15 and 0.85 seems stronger than 
for the lower mass ranges, where it is substantially constant over 
the same redshift range. For the high mass range, the increase in 
the bias going back in time corresponds to a factor ~ 1 .2 (from 
b = 1.40 + 0.03 lob^ 1.62 + 0.18) 

This kind of behavior is generally expected for the clusterin g 
evolution of dark- matter halos of different mass. ( Mo & Whitd 
ll996l:ISheth et al.L i2001; Mo & White, ,2002) . More massive ha- 
los form near the peaks of the density field, whose distribution is 
strongly biased with respect to the overall mass distribution and 
displays a higher clustering amplitude (Kaiser, 1984). Our result 
therefore suggests a direct relationship (as one would naively ex- 
pect), between the total mass in stars accumulated today (or at 
z ~ 1 ) and that of the dark-matter halo within which the bulk 
of the baryonic material was originally assembled. Clearly, to 
be still observed at these redshifts, this proportionality needs 
to have been preserved also through the merging events experi- 
enced by the galaxy. These produce a merging of the dark-matter 
halo and a parallel full coalescence of the baryonic components, 
with the new baryonic remnant still sitting in the center of its 
new, more massive halo. This is not the case when the galaxy is 
accreted by a much larger halo (e.g. a cluster), becoming a satel- 
lite whose own halo will be significantly disturbed and modified 
in mass by dynamical frictio n. This scenario fits well with the 
results of recent simulations dConrov et al.L 120061: IWang et all 
12006 ). that are able to reproduce the observed clustering prop- 
erties of galaxies as a function of luminosity at very different 
redshifts (up to z ~ 3) by simply assuming that each dark-matter 
sub-halo would generate a galaxy characterized by a B-band lu- 
minosity proportional to its original "infall" mass (i.e. mass be- 
fore being last accreted into a larger halo). 

7. Summary and conclusions 

We have used the VVDS-Deep data dLe Fevre et al.L l2005al) to 
perform a first exploration of the dependence of galaxy cluster- 



ing on stellar mass at a redshift approaching unity. This has only 
now become possible at these epochs, thanks to the large number 
of redshifts and extended wavelength coverage available from 
VVDS-Deep that allow us to compute reliable stellar masses. 
We started with a sample of 3218 galaxies with M > 10^ Mq 
in the redshift range z = [0.5 - 1 .2]. We thoroughly investigated 
the completeness limit in mass of the redshift catalogue, using 
mock samples built from the Millennium simulation coupled to 
semi-analytical recipes. We have found that there is a significant 
mass incompleteness induced by the flux-limit of the survey be- 
low M - 10'*' Mq, whose effects on the measured clustering 
become particularly severe below M = lO''-^ Mq. We used the 
mock surveys to quantify this systematic effect and to estimate 
realistic statistical error bars that possibly include "cosmic" vari- 
ance from fluctuations on scales larger than the explored volume. 
With these limitations in mind, we have obtained a series of re- 
sults, that can be summarized as follows: 

- We have measured for the first time a clear stellar mass- 
dependent clustering of galaxies with respect to their stellar 
mass at a redshift as high as z ~ 0.85, with the more mas- 
sive objects being more clustered than less massive ones. In 
particular, the most massive objects show an increase in both 
amplitude and slope of the spatial correlation function. 

- These clustering properties are very well reproduced (within 
1-cr) in the same redshift and mass ranges, by the mock 
samples built from the Millennium run. Together with their 
realistic redshift distribution, this represents a remarkable 
achievement of the models. 

- Comparison of our measurements to the similar ones ob- 
tained (wi th much better statistics) from the SDSS data at 
z ~ 0.15 (iLi et al.L 12006) show evidence for a more signif- 
icant change in the apparent clustering amplitude for low- 
mass galaxies (log(M/M0) < 10) than in the case of the most 
massive objects (log(M/M0) ~ 1 1). The correlation function 
of the latter, in fact, remains roughly constant with time. 

- Assuming a standard cosmology, we computed the expected 
evolution of the root-mean-squared fluctuation in the mass, 
erg, between the two epochs sampled by VVDS and SDSS. 
Using this value, we translated the observed evolution of 
Wp{rp) into the corresponding linear bias, showing that the 
most massive galaxies display an evolution of the bias fac- 
tor, from b = 1.40 + 0.03 at z ~ 0.15 to = 1.62 + 0.18 
at z ~ 0.85. This finding is not unexpected in a hierarchi- 
cal scenario in which the most massive peak s of the mass 
density field collapse earlier and evolve faster dMo & Whitd 
1996). This interpretation supports a scenario in which the 
stellar mass of a galaxy is essentially proportional to the 
mass of the dark-matter halo in which it wa s last the cen- 
tral object, consistent with recent simulations dConrov et al.L 
■2006: Wang et al.,.2006,) . 

The combination of these measurements of clustering as a 
function of ste llar mass and as a function of luminosity at the 
same redshift dPoUo et al.L |2006), together with the analogous 
SDSS results at z ~ 0.15, can provide very important comple- 
mentary constraints to models des cribing how galax ies are dis- 
tributed within dark matter halos dZheng et al.L 12007 ). This can 
help us to understand the evolution of galaxies of different mass 
below z ~ 1 and in particular on the role of mergers within this 
redshift range. 
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